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Abstract. The effect of rotation on the Frequency Ratio Method jMova et alj2005h is examined. Its applicability to observed 
frequencies of rotating y Doradus stars is discussed taking into account the following aspects: the use of a perturbative approach 
to compute adiabatic oscillation frequencies; the effect of rotation on the observational Brunt-Vaisala integral determination 
^ ' and finally, the problem of disentangling multiplet-like structures from frequency patterns due to the period spacing expected 

04 ^ for high-order gravity modes in asymptotic regime. This analysis reveals that the FRM produces reliable results for objects 

OO . with rotational velocities up to 70 km s~', for whic h the FRM intrinsic error increases one order of magnitude with respect to 

' the typical FRM errors given in lMova et alj (2005^. Our computations suggest that, given the spherical degree £ identification, 

the FRM may be discriminating for m = modes, in the sense that the method avoids any misinterpretation induced by the 
presence of rotationally split multiplet-like structures, which reinforces the robustness of the method. However, if £ is unknown, 
such discrimination is not ensured. In order to check the FRM in presence of slow-moderate rotation, we have applied it to the 
three observed frequencies of the slowly rotating (vsin; = 29 km s" ' ) y Doradus star HD 48501. 

Key words, stars: rotation - stars: variables: general stars: oscillations (including pulsations) - stars: fundamental parameters 
Q . - stars: evolution - stars: individual: HD 48501 

■ 1 . Introduction method was applied to the y Doradus star HD 12901, for which 

T ITT n J ^ , a significant reduction of possible representative models of the 

In IMova et al.l 92005) (hereafter paper I), we presented , ■ , 

^ ' r — ; — "r ... . . . ^ . ^ star was achieved. 

>^ . a method for obtaining asteroseismic information of r„, , , , , ^ 

_ , , , , ■„ ■ ^ The FRM, based on the first-order asymptotic ?-mode 

?H y Doradus stars when at least three oscillation frequen- , n 17777^ , „ , . 

. , J T,, . J r 11 J T-niv^ expression given by ITassoull (1198 0) (see Eq. 1 in paper I), 

. . cies are observed. This method, from now on called FRM . , ., f , ■ %■ ,. , . . , 

,„ r. • . J 1 1^ • -11 ^ 1 ^ • ■ IS built under the assumptions of adiabaticity and non ro- 

(Frequency Ratio Method), IS particularly useful for improving . . ^ 

, , , , 1 • 11 1 • ^ J tation. However, recent spectroscopic measurements of 59 

our knowledge on the pulsational behaviour of y Doradus stars ^ , ,. , l , |. — . .1 , 

, 1 1 11 • n , ■ A n ■ r , 7 Doradus candidates C Mathias et alJ2004D show rotational ve- 

through the following aspects: first, the identification of the ... ^ . — , _i 

J, p r , J r ■ loci ties (vsin; measurements) varying from 10 to loOkms . 

radial order n and the mode degree t of observed frequencies ^, , , , • , , ^ ^ ^ ^ , 

, , , -1 -1 .c • • J 1 This corresponds to rotational periods around 0.4 to 4-5 days, 

and second, the possibility of constraining models from an . , . ,• ,, 

J J. . , „ . a / I • 1 1 1 in the limit of applicability of a perturbative analysis used here 

additional stellar quantity, J (the integral of the buoyancy ,. , . .,, . ^ . 

T, , / • 1 J ^1 ^ 11 J- to compute adiabatic oscillation frequencies. 

or Brunt-Vaisala frequency weighted over the stellar radius ^ i ^— ' — i j 



, A- \ c u * • t u ■ » ^ A It has been shown by [Michel et al. ( 1 19991) and references 

along the radiative zone). Such constraints can be interpreted ' ' ^ 

, ■ 1 u ui therein, that photometric parameters are sensitive to rotation. In 

as an increase of the classical observables, ^ ^ 

particular, it modifies the determination of the position of the 

y - y(Z L T (f ■ 3) ^'■^^ ^'^ '■^^ diagram. Since the FRM results are sensitive to 

' foj ' ' such variations, the rotation must thus be included in the set of 

working parameters 



where Z represents the relative metal abundance; g the surface 
gravity, L the luminosity; T^s the effective temperature, and 



finally, /o,, and f^j the observed frequencies. In paper I, the ^ - ^n{Z, g, L,Tsf[; —^,3,D.), 

'■ '■ Jo.j 
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* Associated researcher at (2) where Q represents the rotational velocity of the star. 
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Concerning the theoretical adiabatic oscillations, special 
care must be taken when the oscillation frequencies are of 
the order of magnitude of the rotational frequency of the 
star, i.e. when ctq ~ cr/3. For such cases the perturba- 
tive method for computing adiabatic frequencies may start 
to fail, and other non-perturbative theories should be used. 
In t his c ontext, several works can be found in the litera- 
turejDintrans & RieutordI (I2OO0I) . iRieu tord & Dintrans (2002), 
IPintrans & BrandenburgI ll2004l) . in which non-radial gravity 
modes of a typical y Doradus star are studied by means of poly- 
tropic models. 

In addition, the effects of rotation on the oscillation spectra 
may introduce additional uncertainties to the FRM. The prob- 
lem arises when the classical g-mode frequency pattern may 
be misidentified with multiplet-like splitting rotationally in- 
duced structures. Ideally, when considering complete asymp- 
totic spectra of y Doradus stars, in a first approximation, 
the asymptotic g-mode pattern, whose periods are regularly 
spaced, can be easily disentangled from multiplet-like struc- 
tures, which are regularly spaced in frequency. However in 
practice, unfortunately only a few oscillation modes are ob- 
served, and therefore, such direct discrimination is not accom- 
plishable. 

The paper is organised as follows: In Sect. |2] the equilib- 
rium models are described. The domain of validity of the FRM 
when applied to rotating y Doradus stars is analysed in Sect.|3] 
In Sect.|3 the problem of multiplet-like structures in the frame- 
work of the FRM is discussed. An application of the FRM 
for the rotating y Doradus star HD 48501 is given in Sect. |5] 
Finally, conclusions are summarised in Sect.|6l 

2. The models 

For the purpose of constructing representative models of 
y Doradus stars the evolutionary code CESAM (Morel 1997) 
is particularly adapted. Models are built with a precision op- 
timised to compute oscillations, i.e. around 2000 mesh points 
for the equilibrium model mesh grid, given in the basis of B- 
splines. To take into account a first order effect of the rotation, 
equilibrium equations are modif ied in the manner described 
in Ucippenhahn & WeigerJ lll990l) . In particular, the spherical 
symmetric contribution of the centrifugal acceleration is in- 
cluded by means of an effective gravity 

gss^ g- -Acir) , 

where g represents the local gravity component, r the radial 
distance and, 

the centrifugal acceleration of matter elements at a distance r 
from the centre of the star. This spherically symmetric contri- 
bution of the rotation does not change the shape of the hydro- 
static equilibrium equation. Although, the non-spheric compo- 
nents of the centrifugal acceleration are not considered, they 
are included as a perturbation in the oscillation computation. 
The total angular momentum of models is assumed to be glob- 
ally conserved along the evolution of the star. 



The physical parametrisation has been adapted for inter- 
mediate mass stars (see paper I). In particular, all equilib- 
rium models used in the present work are computed with : 
cml - Im/Hp = 1-8, for the convective efficiency and dgv - 
l„v/iip = 0.2 for the mixed core overshooting. The Hp corre- 
sponds to the local pressure scale-height; /„, and dov represent 
the mixing length and the inertial penetration distance of con- 
vective bulbs respectively. 

3. The error determination 

The study of the influence of rotation on the FRM constitutes 
a rather complex task. We aim at identifying possible uncer- 
tainties coming from rotation effects and their relevance for the 
FRM. Such analysis is here focused in three important aspects: 
first, the regime of application of the perturbative method for 
obtaining theoretical adiabatic frequencies. Second, the quan- 
tification of the error committed when considering the classi- 
cal analytical expression for asymptotic g modes to study fre- 
quency ratios of rotating y Doradus stars. And third, the influ- 
ence of rotation on the Brunt- Vaisala integrals 3. 

3.1. The perturbative approach for adiabatic oscillation 
frequencies 

The limit for using the perturbative approach to compute adia- 
batic oscillation frequencies is given by the relation cr/crn » 1, 
i.e. when the rotational frequency of the star (crn) is negli- 
gible with respect to the the oscillation frequency of modes 
(cr). In the asymptotic regime, g modes are characterised 
by their low frequency (high radial order n). According to 
Ipintrans & Rieutord ( 2000), the most perturbed modes by rota- 
tion are likely those with cr < 2 ctq, where 2 ctq corresponds to 
the Coriolis frequency. Considering the typical rotational peri- 
ods of y Doradus stars, low frequency gravity modes in asymp- 
totic regime are strongly affected by rotation. 

The present analysis is thus restricted to slow and moder- 
ately rotating y Doradus stars, that is, those with rotational ve- 
locities smaller than 50-70 km s"'. Considering representative 
models with masses in the range of 1 .4 - 1 .7M0, the limits of 
the perturbative approach can be established to have «=[30, 40] 
for modes with ( < 2, expected to have significant amplitude 
not smeared out by cancellation effects. 

3.2. The analytic expression for the asymptotic g 
mode frequencies 

It is convenient to remind that the analytic expression for the 
asymptotic g modes used in the FRM is obtained under the 
hypothesis of non-rotation (see paper I). We attempt here to 
determine the accuracy of this expression when rotation effects 
are considered. In this framework, it is conceivable to define an 
equivalent form of the analytic expression for the asymptotic g 
modes as: 

where the analytic form of integrals remains identical to 3 
(see paper I), and whose values are assumed to be different to 
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Fig. 1. Relative eirors of the asymptotic frequencies as compared to the theoretical predictions for a complete oscillation spectrum 
obtained from a model with Q ~ 50km s '. Left panel represent relative errors obtained from oscillations frequency ratios 
coiTected for the effect of rotation. Right panel shows errors without correcting for the effect of rotation. 



those obtained without considering rotation effects. The f(n, €) 
functions are a priori unknown and assumed to be only depen- 
dent on { and n numbers. Similarly as done in paper I, two given 
g modes, ctq,, q and cr„, q, are considered, where a\ and ai rep- 
resent the pairs (ni, ^i) and {ni, ^2) respectively. Assuming the 
same mode degree t\ - I2 - (, we make the hypothesis that 
their eigenfrequencies are approximated by Eq. If the star 
rotates slowly (cf. Sect. 13. U . f{n, {) can be assumed to be: 



fin. 



and Eq.[2can thus be rewritten as follows: 



(2) 



(3) 



(n+l/2)7rj^ 

In this way, the analytic form proposed bv lTassoull ill 9801) for 
non-rotating stars is kept and consequently, the formal structure 
of the FRM remains unaltered. This allows us to relate crcn,si 
and ctq.^ n through the following ratio 



O-Q. 



«2 + l/2 \ 



(4) 



To obtain this expression, we proceed similarly as in paper I 
(Sect. 3). It can be shown that the small differences in the outer 
turning point location can be considered as negligible as far 
as the calculation of the Brunt-Vaisala frequency integral is 
concerned. 

Provided theoretical frequencies and the corresponding n's 
ratio estimate for a given model, the accuracy of Eq0]can be 
defined as 

'(r,\ (f{njj)\ 



(5) 



■J /a \ /(«<■' ^)/n 
Figure [1] (left panel) shows such errors obtained from a theo- 
retical oscillation spectra of a 1 .40 model. It can be noticed 
that errors reach 0.013 for ^ = 1 modes, and 0.021 for ^ = 2 



modes, for rotational velocities around 50-60 km s ' . These re- 
sults represent around one order of magnitude larger than the 
e/j o errors (right panel), obtained without taking into account 
corrections for rotation in frequencies. As a consequence, the 
quantity of valid sets of natural numbers (n,) increases (see pa- 
per I, Sect. 3). Such increase can be of the order of 10 addi- 
tional sets for { - 1 modes, and even 100 sets for £ - 2 modes. 
Nevertheless, as explained in paper I each additional possible 
set, implies an extra constraint for the Brunt-Vaisala integrals 
which must be fulfilled. 

3.3. The effects of rotation on the 3 integrals 

This study must be completed with the analysis of the be- 
haviour of 3 integrals when rotation is taken into account. 

As explained in paper I, in the FRM framework the Brunt- 
Vaisala frequency integrals are calculated in two different man- 
ners: from equilibrium models and directly from frequency ra- 
tios as deduced from the analytic expression for asymptotic g 
modes (Eq. 1^. The characteristics of pseudo-rotating mod- 
els (cf. Sect. O allow us to consider equivalent internal struc- 
ture for rotating y Doradus stars only modified by an effective 
gravity. Therefore, the Brunt-Vaisala frequency integrals can 
be calculated as for non-rotating models, that is 



cr 



n,t,n 



'th.n = 



fin J) 



(6) 



These integrals are dependent mainly on the mass and metal- 
licity of models (Fig. 5 in paper I). For rotating 7 Doradus stars 
it is plausible to consider also a Q dependence 



= Jth(M,Z,Q). 



(7) 



For the purpose of investigating the influence of rotation on 
Brunt-Vaisala integrals, we consider 1.4 M© pseudo-rotating 
models as described in Sect. |2 Models are evolved from the 
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Fig. 2. Variation of theoretical Brunt-Vaisala frequency inte- 
gral as a function of the evolutionary stage (represented by the 
central Hydrogen fraction Xc) for a I.4M0 model computed 
with solar metallicity. Symbols represent models with different 
initial (ZAMS) rotational velocity. 



Table 1. Characteristics of selected models in Fig.|3j present- 
ing a common 3th ~ 700 //Hz, and evolved from the initial (at 
ZAMS) rotational velocities 25,50,100 and 150 kms^^ (Qi, 
f22, ^3 and ^\ respectively). The following stellar parameters 
are given: Teff, the effective temperature (in logarithmic scale); 
g, the surface gravity (in a logarithmic scale); p, the mean stel- 
lar density (in g cm^); R/Rq, the stellar radius (in solar units); 
Xc, the central Hydrogen fraction; Q, the rotational velocity 
(in kms '), and finally, 3ih, the theoretical Brunt-Vaisala fre- 
quency integral (in juHz). 





n4 




n2 


Hi 


log Teff 


3.81 


3.81 


3.82 


3.82 


log g 


3.96 


3.96 


3.97 


3.97 


P 


3.79 


3.75 


3.88 


3.84 


R/Rq 


2.18 


2.19 


2.17 


2.17 


X, 


0.32 


0.30 


0.29 


0.28 


n 


132.50 


87.87 


42.18 


20.99 


3th 


697.81 


703.31 


698.15 


702.08 
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Fig. 3. Theoretical Brunt-Vaisala integral as a function of the 
rotational velocity Q. Filled grey circles represent the same 
1 .4 Mq models evolved with the same initial rotational veloci- 
ties considered in Fig.|2] Filled black circles represent selected 
models with the 3th values closest to 700 juHz (see details in the 
text). 



ZAMS to the subgiant region with four different initial rota- 
tional velocities: 25, 50, 100 and 150 kms"'. For these mod- 
els, the evolution of 3th,n is depicted in Fig. |2] as a function 
of the core hydrogen fraction X^. Small variations (< 10/vHz) 
are observed. As expected, the effect of considering an effec- 
tive gravity does not modify the behaviour of the theoretical 
Brunt-Vaisala integrals along the stellar evolution. 

In practice, initial rotational velocities are unknown. 
Therefore, it is worth to analyse the characteristics of mod- 
els presenting a similar Brunt-Vaisala integral but having been 
evolved with different initial rotational velocities. In Fig. |3] 



the evolution of 3th,n integrals is shown as a function of the 
rotational velocity. Each track corresponds to evolved models 
(from left to right) with Ozams =150, 100, 50 and 25 kms"'. 
From each set of models, those with 3th,n~ 700 //Hz (black cir- 
cles) have been arbitrarily selected. The analysis of the cor- 
responding stellar parameters given in Table ^ reveals simi- 
lar characteristics for the four models. Small differences are 
mainly due to the different evolutionary stage of models. In 
particular, the +5 /iHz dispersion in 3th,n values are due to 
slight differences in their integrands (Fig.|4}. As expected, ro- 
tation does not discriminate models with respect to their 3th,n- 
Nevertheless this ensures the coherence of the FRM, since rota- 
tion does not introduce additional uncertainties from the point 
of view of the 3th behaviour 

Le us consider now the observed Brunt-Vaisala frequency 
integrals. The analysis of uncertainties of such integrals is not 
trivial. According to Eq.[2 these integrals can be expressed as 



O" n,C,Q. 

fin J) 



(8) 



From this expression and Eq.|6l it is possible to define the error 

associated to e^^.n as 



|o-«/ ,n - f(n, Jth,nl 
o" n,e,n 



(9) 



Respectively, the error associated to the Brunt-Vaisala fre- 
quency integral can be written as 



th,n ■ 



Jobs.nl 



hh,n 



(10) 



Multiplying both sides of Eq.^Jby /(n, O' , and then consider- 
ing the Eqs.|S]and Eq.|5] both errors can be related as follows: 



This operation is licit since f(n, ^) > by construction. 
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Fig. 4. Brunt- Vaisala frequency weighted over the stellar radius (in solar units) as a function of the radial distance r (normalised 
to the stellar radius R), for the model with the rotational velocity Qj (see Table Right panel shows the relative difference 
between the Brunt-Vaisala frequency weighted over the stellar radius of models with Q.\ and ^4 respectively (see Table[l)- 



£a-,Q. 



/(«, €) Jth.n 



(11) 



According to Eq. |5] e^r.n can be discussed as a function of the 
sign of 



O" n,t,Q. 



Jth.n 



+ea-,Q. if ^^-^ 3th,n > 1 
cr n,e,si 

-fo-.n ir Jth.n < t 



(12) 



where the trivial case of fo-n = is obviously omitted. The 
right hand side of Eq.^Jcan thus be written as 



/(«, f) Jth,n 



= 1 + e„,f,n , 



(13) 



This allows us to quantify the error in Uobs.n integrals as a func- 
tion of fer n by replacing the left hand side of Eq. [^in Eq. [01 
yielding 

ej.n = , (14) 

1 ± e<T,n 

This means that ej n errors can thus be considered of the same 
order as ea-si- Such result was somehow anticipated in pa- 
per I from numerical computations. Here, although the relation 
given by Eq.^^has been obtained for cr„/_n and Jobs.n errors, 
it can be easily shown that is generally valid whatever consid- 
ered errors and /(n, €) values. Similarly as it happens for non- 
rotating models, 63 n ~ 10"^' which implies an error of 10"^ in 



the FRM predictions may fail. In particular, rotationally split 
multiplets (or any of their m + ^ components) may be con- 
fused with the large separation pattern expected for high-order 
gravity modes in asymptotic regime. 

In order to investigate this, an observed oscillation spec- 
trum of a given rotating y Doradus star has been simulated. 
This simulation has been performed using a pseudo-rotating 
model (as described in Sect. |2ji with a mass of 1 .70 M©, and a 
rotational velocity of Qq = 26.8 kms '. Computed with solar 
metallicity, it is located in the main sequence (age 1.3 Gyr). It 
presents an effective temperature and a surface gravity (log 
cgs) of 6918 K and 3.94 respectively. Its theoretical Brunt- 
Vaisala frequency is Jth,n= 667.8 |/Hz. 

Adiabatic eigenfrequencies have been computed 
from this model with the oscillation code FILOU (see 
iTran Minh & LeonI Il995'l) . In this code, effects of rotation 
up to sec ond order on t he oscillati on frequenc i es hav e been 
included d Suared l2Q02l) following ISoufi et alJ lll998l) . Due 
to the nature of the observed frequencies, the theoretical 
oscillation computations are extended to high order g modes 
(-40 < n < 1), limiting our investigations to modes with 
degrees { < 2, since cancellation effects are important for the 
visibility of modes with { >3. 

To cover all the possible situations where multiplet-like 
structures may contaminate the FRM results, the following sce- 
narios have been explored: 



4. The problem of multiplet-like structures 

When rotational splittings are of the same order as the period 
spacing of high radial order g modes, i.e., 

|o"«/,m - cr„/,,„'| ~ |cr„/_o - cr„'/,ol , (15) 



1 . the observed frequencies belong to any rotational multiplet; 

2. they belong to modes with diflFerent radial order n and az- 
imuthal orders m; 

3. they are identified as modes with different (n) but with m - 
0. 
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Table 2. List of possible (n, J) identifications provided by the 
FRM applied to the theoretical frequencies belonging to the 
rotational (« - 22, { - I) triplet (more details in text). The first 
three columns represent the resulting radial order identification. 
Last two, list the corresponding observed 3 when assuming i - 
1 and 2 respectively. 





H2 


/i3 






ZD 


li 


21 


nr\A inc\ 

/04.J /y 


4U6.6/3 


27 


24 


22 


734.353 


423.978 


28 


25 


23 


764.326 


441.284 


36 


32 


29 


974.141 


562.421 


37 


33 


30 


1004.115 


579.726 


38 


34 


31 


1034.089 


597.031 


39 


35 


32 


1064.062 


614.336 


40 


36 


33 


1094.036 


631.642 


42 


37 


34 


1124.009 


648.947 


43 


38 


35 


1153.983 


666.252 


44 


39 


36 


1183.957 


683.557 


45 


40 


37 


1213.930 


700.863 



Table 3. Same as Table |2]but obtained for the case of modes 
with different radial and azimuthal orders (more details in the 
text). 



Ml 


«2 


"3 


?obs(^ - 1) 


3obs(^ = 2) 


20 


16 


13 


494.564 


285.536 


25 


20 


16 


614.458 


354.757 


26 


21 


17 


644.432 


372.063 


31 


25 


20 


764.326 


441.284 


32 


26 


21 


794.300 


458.589 


37 


30 


24 


914.194 


527.810 


38 


31 


25 


944.168 


545.115 


40 


32 


26 


974.141 


562.421 


42 


34 


27 


1034.089 


597.031 


43 


35 


28 


1064.062 


614.336 


44 


36 


29 


1094.036 


631.642 


45 


36 


29 


1094.036 


631.642 


46 


37 


30 


1124.009 


648.947 


47 


38 


31 


1153.983 


666.252 


48 


39 


31 


1183.957 


683.557 


49 


40 


32 


1213.930 


700.863 



4A. Scenario 1 

For the first scenario, let us examine, for instance, the observed 
(n, €) — (22, 1) triplet with frequencies (in //Hz) 

cr_i = 14.756, o-o = 13.493, cr+j = 12.045, 



and with frequency ratios 

^. 0.8927, ^ 



0.9144, — =0.8163. 

o--\ 



For these frequency ratios, the FRM predicts a list of possi- 
ble (n, 3) identifications indicated in Table |2l As discussed in 
Sect. 13.21 this list has been constructed by assuming an error of 
+ 1.3 10"^ for ep. An estimate of 3obs,n integrals is also given for 
the two most relevant mode degree values: I - \ and 2. The ob- 
served set is obviously not predicted since n\ -n^ -n^. When 
analysing the list of possible Jobs.n integrals, predictions ex- 
ceed the error limit (around 2%) when assuming I -\. That is, 
for this known value of the mode degree, none of the potential 
sets would be misinterpreted. The same exercise repeated with 
other I - \ multiplets within the same spectral range yields 
similar results. In contrast, the list of Brunt-Vaisala integrals 
predicted assuming 1-2 contains Jobs.n values close to the 
theoretical one. In particular, the set {n\,n2,ni) - (43,38,35) 
is related to 3obs,n(^ = 2)=666.252/iHz, which represents an 
error of 0.23% approximately with respect to the 3th- This case 
would thus be relevant when no additional information on the 
mode degree is provided. In such a situation, the discrimina- 
tion between a possible multiplet-like structure and asymptotic 
g mode pattern is not ensured. 

4.2. Scenario 2 

The second scenario implies different radial and azimuthal or- 
ders. To check it, the £ - I modes with (n, m)=(24,-Hl), (22,0) 



and (19,-1) have been arbitrarily selected. Their corresponding 
frequencies are 

cri9,_i = 16.727, 0-22,0 = 13.493, o-24,+i = 10.979, 

for which the corresponding frequency ratios are 

o"+i cm cr.i 

— =0.8137, —= 0.8066, —= 0.6563. 

(To o-_i cr_i 

In this case, the application of the FRM to these frequency 
ratios yields the list of possible («, 3) identifications given in 
Table |3] Once again, the solution is not predicted by the FRM. 
Moreover, the 3obs,n predictions for ^ = 1 modes are out of the 
error limit. However, as discussed in Scenario 1, when consid- 
ering ( - 2 modes, the theoretical Brunt-Vaisala integral is 
predicted, at least by one of the sets: («i, n2> «3) = (47, 38, 31). 



4.3. Scenario 3 

The third and last scenario is constituted by three centroid 
modes (m - 0), with different radial order Proceeding as in 
Scenarios 1 and 2, the frequencies of the { = I modes with 
« = 19, 22 and 24 are 

cr,9 = 15,4633, cr22 = 13.4928, cr24 = 12.4232, 

and their frequency ratios are 

— = 0.9207, — = 0.8726, — = 0.8034 . 

(T22 Cri9 (Tig 

The (n, 3) possible identifications given by the FRM applied to 
these frequencies are listed in Tabled] In contrast to the pre- 
vious scenarios, when considering centroid modes the solution 
is found by the FRM. Not only the observed radial order set is 
included in the list of possible solutions (and thereby the mode 
degree {), but also its related Jobs.n value predicts the theoreti- 
cal Brunt-Vaisala integral within the error Umits. 
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Fig. 5. Frequency differences between the components of the observed modes proposed in Scenario 3 (Sec. 14. 3t and those of 
theoretical multiplets as a function of the radial order n. Left and right panels show frequency differences calculated considering 
{ - \ and { - 2 modes respectively. Symbols represent the frequency differences calculated for the three observed I - 1, m = 
modes: circles for crig, triangles for 1x22 and finally, stars for 0-24. Azimuthal orders should be read as going from negative to 
positive (m - +|^|) from left to right respectively. Filled circles represent modes for which bcr 5 0.5 (grey) and 6^ < 0. 1 (black). 
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Table 4. Same as Table|2]for the case of modes with different 
radial order but the same (m - 0) azimuthal order (more details 
in the text). 



n, 


"2 


«3 


Jobs,n(^ - 1) 


Jobs,n(^ - 2) 


24 


22 


19 


674.405 


389.368 


25 


23 


20 


704.379 


406.673 


26 


24 


21 


734.353 


423.978 


34 


31 


27 


944.168 


545.115 


35 


32 


28 


974.141 


562.421 


36 


33 


29 


1004.115 


579.726 


37 


34 


30 


1034.089 


597.031 


39 


36 


31 


1094.036 


631.642 


40 


37 


32 


1124.009 


648.947 



0.060 F 



0.050 - * 



0.040 r 

N : 

^ O.OSO r ★ 
e * 

I ; • , 

^ 0.020 r 



0.01 O'r 



0.000 F, , , ,^ I , , ,. 1^ » I ^, , , I F 

27. to 27.00 26.90 26. BO 26.70 26.60 
f! (kra/s) 



4.4. Considering different rotational velocities 

The previous discussion considers only one rotational veloc- 
ity. However, usually, it cannot be determined from observa- 
tions neither the specific surface rotational velocity nor the an- 
gle of inclination of the star but the value of vsin/. Since fre- 
quency patterns are Q. dependent, it is plausible to imagine a 
given rotational velocity for which Eq.[^ holds. To investi- 
gate this possibility, 7 models with different rotational veloc- 
ities have been considered: three models with Q < Qq, three 
models with Q > Qq and finally the model used for the the- 
oretical test which rotates at Qq = 26.8 km s"'. The range of 
rotational velocities has been chosen as to produce models rep- 
resentative of a typical y Doradus star The HR diagram error 
boxes have been obtained from typical multicolour photometry 
errors of y Doradus stars. 

We are looking for any trend indicating that the observed 
frequencies crig, 0-22, and 0-24 may be misidentified with a ro- 
tational triplet or with components of any { = 1,2 multiple ts. 
To do so, we have calculated the frequency difference da- be- 
tween the observed modes and those of theoretical multiplets 
for the 7 models considered here. In Fig.|5] such differences are 
shown for the reference model (central panels), the fastest one 
(top panels) and the slowest one (bottom panel). In the case of 
{ - I triplets (left panels), the smallest \Sa-\ values vary from 
0.1 (black filled circles) to 0.5 (grey filled circles) approxi- 
mately, for a set of 15-20 modes. As expected, for £ = 2 modes 
(right panels) da- decreases, and the number of modes with fre- 
quencies close to the observed ones increases. In this case, the 
smallest |<5o-| values vary from 0.001 fiHz (black filled circles) 
to 0.3 /iHz (grey filled circles) approximately, for a set of 35-40 
modes. Discarding the exact solution for crig, a-22, and 0-24 for 
which 6cr - 0, these frequency differences exceed in at least 
three orders of magnitude {{ - 2 results)^ the observational 
frequency error typical for y Doradus stars. 

In order to consider the smallest differences for the 7 
seven models, the minimum da- value found for each model 
is depicted in Fig. |6l as a function of the rotational velocity. 
Although the mode degree is not here specifically identified, it 

^ For £ = \, 6a- differences reach five orders of magnitude that of 
the observational frequency error. 



Fig. 6. Variation of the minimum frequency differences consid- 
ered in Fig. |5] with the rotational velocity. Same symbol repre- 
sentation as in Fig.|5]is adopted. 



is found that most of selected modes are { -2, result somehow 
expected since, according to Fig. |5] the smallest differences 
corresponds to i -2. Obviously, the exact solution {6a-,mm - 
at n = Qo) correspond to ^ = 1 modes. As can be seen, the 
5cr,mm distribution does not show any specific trend with the 
rotational velocity and, quantitatively, the results obtained are 
similar to those shown in Fig.|5l This result does not exclude 
the possibility of misidentification, however it shows it is un- 
likely for the range of models representative of the observed 
star. 

This theoretical test shows that, provided the mode de- 
gree { is known, the FRM may be discriminating with respect 
to the multiplet-like structures. In such cases, coherent FRM 
solutions are restricted to centroid mode {m - 0) frequen- 
cies^. When no additional information on i is given, the FRM 
does not discriminate multiplet-like structures. Nevertheless, as 
shown in this exercise, such a situation is unlikely. 

5. A particular case: the y Doradus star HD 48501 

The y Doradus star HD 48501 is particularly suitable to il- 
lustrate the previous theoretical exercise applied to a real 
star. Having a vsin/ ~ 30kms"\ it can be considered as a 
slow rotator, fulfilling thus one of the conditions discussed in 
Sect.lSl Ac cording to a recent frequency analysis reported by 
lAerts et aLi (2004) the star pulsates in at least three modes. In 
this paper, a discussion has been opened on whether these fre- 
quencies could belong to a f = 1 triplet-like structure or to the 
period spacing expected for high-order gravity modes. 

5.1. Fundamental parameters 

The stellar fundamental parameters o f HD 48501 have been 
determined by applying TempLogG ( lKupka&Brunttll200ll) 

^ Such result is obviously restricted to stars fulfilling the conditions 
and errors discussed in Sect.|3| 
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Table 5. Fundamental parameters of HD 48501 taken from 
the literature. From left to right: effective temperature (in 
K); surface gravity (cgs, in a logarithmic scale); metallic- 
ity and finally the projected rotational velocity (in kms 
Superindex indicat es data were ob t ained from the follow- 
ing re ferenc es: (a ) Kupka&Bruntt (200lk; (^) [Ducret et al 
1 20021) : (c) ICavr el de Strobel e t aL ^199% : (t/)^ 
1 2004 and finally (e) Rover et alfi Eooa ^ 



Table 6. Photometric data for the y Doradus star HD 48501. 
Observed frequencies, in bo th cycles per day and juHz, are 
taken from lAerts et alJ ll2004 . 



^erts et al 





log g 


[Fe/H] 


vsin; 


(K) 


(dex) 


(dex) 


(kms-') 


6984" 


3.92" 


-0.12" 


29'' 


7079* 


4.49* 


-0.10* 


47<- 


7099'' 


3.96'' 


+o.or 





to the Stromgren-Cra wford photometry hsted in the Hauck- 
Mermilliod catalogue jHanck & Memiillindll 998l) . In this cat- 
alogue, no y6 index value is given for HD 48501, which was 
instead obtained from Handler (1999). The code classifies this 
object as a main sequence star in the spectral range region FO- 
G2. The resulting physical parameters are listed in Table|5l 

bupret et al. ( ,2002.) gives also physical par ameters for 
HD 48 501 based on the Geneva ph otometrv of lAerts et alJ 
(12004') and the calibrations given bv lKunzliet all ( 11997 ') for 
the Geneva photometry of B to G stars. In his study, theo- 
retical stellar models for these stars do not take into account 
their surface gravity determinations, for being too high, cor- 
responding to models below the ZAMS. Physical parameters 
were also retrieved from the fifth edition of the catalogu e 
of [Fe/H] determinations by 'Cavrel de Strob el et al ] (Il997h . 
which includes [Fe/H] and atmospheric parameters (reff, log g) 
obtained from high resolution spectroscopic observations and 
detailed analysis, most of them carried out with the help of 
model-atmospheres. The aforementioned physical parameters 
are listed in Table |5l 

Values for the vsin/ of were found in 'Rov er et alJ tOOf). 
computed from spectra c ollected at O bservatoire de Haute- 
Provence (OHP) and by lAerts et all (|2p04), derived from 
a cross-correlation function analysis of spectroscopic mea- 
surements wi t h the CORALIE spectrograph. Very recently, 
iMathias et all (l2004 have pubUshed the results of a two-year 
long high-resolution spectroscopy campaign, monitoring 59 
y Doradus candidates. In this campaign, more than 60% of the 
stars presented line profile variations which can be interpreted 
as due to pulsation. 

5.2. Locating HD 48501 in tlie HR diagram 

As usual, the modelling of any star is conditioned by its stel- 
lar parameters, delimiting its position in the HR diagram with 
the corresponding photometric errors. In the present case, the 
photometric error box (see Fig.0 has been delimited from fun- 
damental parameters (Table |5} obtained from the Stromgren- 
Crawford photometry (boldface). For these values, a standard 
accuracy of 0.2 dex for the surface gravity and 200 K in T^s has 





(c/d) 


0"Hz) 


fl 


1.095 


12.663 


fu 


1.199 


13.880 


/ill 


1.290 


14.937 



been considered. In addition, a typical correction for the effect 
of rotation is included as well. For main sequence intermedi- 
ate mass stars, this represents +200K approximately (details in 
[Michel et al. 1999; Perez Hernandez et al. 1999), and more re- 
cently in lSuarez et alJ(l2002l) . In the same figure, several evolu- 
tionary tracks are drawn to illustrate the wide range of different 
models that can be considered as representative of HD 48501. 
In particular, the location in the HR diagram of evolutionary 
tracks can be modified significantly using different rotational 
velocities. This becomes critical when different values for the 
rotational velocity and the metallicity are taken into account. At 
present, there is no additional information constraining the an- 
gle of inclination of HD 4850 1 . Therefore, a wide range of rota- 
tional velocities should be considered. However, in the present 
case, in order to respect the hypothesis of applicability of the 
FRM, we restrict our modelling to low rotational velocities. In 
Fig El the upper continuous line and the lower dashed line rep- 
resent, respectively, the evolutionary tracks of a 1 .76 M© non 
rotating star model and a 1 .65 M© pseudo-rotating star model. 
These tracks illustrate somehow the limits for the mass and ro- 
tational velocities predicted by the modelling. 

Four diff'erent masses have been considered: from bottom 
to top, 1.65Mo {D models), 1.70 Mq {B models) and 1.75Mo 
(A models) computed with solar metallicity. To cover the ob- 
served vsin/ found in the literature (see Table |5}, the follow- 
ing rotational velocities have been used: 30 kms"' {A[, Bi) and 
50 kms"' (A2, B2), corresponding to an angle of inclination of 
90° and 37° respectively. 

5.3. Applying the FRM to tine observed frequencies 

The three observ ed frequenc i es (T able |6t of HD 48501 have 
been taken from lAerts et al] (12004) . The corresponding fre- 
quency ratios are 

^ = 0.9123, ^ = 0.9292, ^ = 0.8478 . 
/ii /ill /ill 

All possible natural number ratios up to n = 60 are then calcu- 
lated. As described in Sect.|3 we assume an error of + 1.3 ■ 10"^ 
for the calculation of possible n sets, instead of the standard er- 
ror of +5 ■ 10"^ used under the assumption of non-rotation (cf. 
paper I). The results are detailed in Table for which only 
dipoles ({ - 1) and quadrupoles ({ - 2) have been considered 
up to n < 50. 

Within the aforementioned errors, seven possible n sets are 
obtained. Although, most of them (5 of 7) are predicted as ^ = 2 
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5. .90 3.88 3.86 3.84 3.82 3.80 3.78 



Fig. 7. HR diagram containing evolutionary tracks and models loci considered in this work. The continuous line represents 
the evolutionary track of a 1.76 non rotating star model. Respectively, the dashed line represents an evolutionary track of 
an 1.65 Mo pseudo-rotating star model, considering a rotational velocity of 50kms '. A, B and D correspond to the location 
of 1.75, Mq, 1.70, Mq and 1.65 pseudo-rotating models respectively (a solar metallicity is considered). Subindexes 1 and 
2 represent values of rotation velocities at the stellar surface of 30kms"' and 50kms"' respectively. C represents models of 
1.50Mo, with a metallicity of Z = 0.07 and a rotational velocity of 50kms Models are displayed with filled circles. Finally, 
the error box considered for HD 48501 is also depicted. 



modes, the and fe are predicted as ^ = 1 modes. This latter 
possibility would thus support the results obtained with mul- 
ticolour photometry given by 'Aerts et aP ( l2004l) . The corre- 
sponding Jobs.n integrals (see Tabled are compared with theo- 
retical predictions in the Jth-7^eff diagram given in Fig.|8l When 
analysing such comparison for ( - \ possibilities (shaded 
boxes), no successful results are afforded, i.e. the models ful- 
filling the («,, U) predictions are not consistent with the observa- 
tional constraints. In particular, the observed metallicity and the 
coiTesponding 3th,n predictions are in clear dissonance. In par- 
ticular, the FRM would predict models with a metallicity in the 
range of [Fe/H]=[-0.4,-0.6]. In contrast, when considering pos- 
sible sets with { - 2, the FRM predicts (according to the Jobs.n 
integrals) models with the observed metallicity [Fe/H]> 0.1, 
in particular for t\, t2 and tj, sets. Such solutions would corre- 
spond to models with masses from 1 .4 to 1 .6 Mq in the main 
sequence. 

5.4. The observed frequencies as belonging to a 
rotational split I = 1 triplet 

The previous results suggest, if ^ = 1 confirmed, the alternative 
of considering the observed frequencies as belonging to one or 
more multiplet-like structure split by rotation. 



Table 7. List of selected natural numbers sets associated to the 
observed frequency ratios of HD 48501. For each set, the cor- 
responding observed Jobs.n is given in juHz. Finally, an estimate 
of the spherical order I is also given in column 4. 
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In Fig. |9] the theoretical evolution of adiabatic frequen- 
cies is given for two representative models of HD 48501. The 
region where the observed frequencies (continuous horizon- 
tal lines) can be identified is around the triplet g22, that is a 
n - 22, g mode. This figure shows that theoretical predictions 
fit better for higher masses (right panel) and for evolved mod- 
els (lower effective temperatures). Two rotational velocities 
are considered: 30 km s" ' , and 50 km s ' , represented by dotted 
and dashed lines respectively. For the low-mass model (Fig.|5J 
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Fig. 8. Theoretical Brunt-Vaisala frequency integral as a function of the effective temperature for representative models of 
HD 48501 with masses in the range of M = 1.2-1.4M0, for [Fe/H] = -0.4 and -0.6, and M = I.2-I.6M0 for [Fe/H] = 0.0, -0.4 
and -0.6. Shaded areas represent uncertainty (3obs,7^eff) boxes of sets ti-t^ and 14, fe as given in Table0 
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Fig. 9. Evolution of adiabatic oscillation frequencies for the 1 .65 M© and 1 .70 Mq models within the HD 48501 photometric error 
box (left and right panels respectively). Only the three components of the g2i triplet are displayed. Dotted lines represent the 
frequencies obtained from models with a rotational velocity of 30kms and dashed lines correspond to those obtained from 
models with a rotational velocity of 50 km s ' . Finally, the observed frequencies are represented by continuous horizontal lines. 



left panel), the predicted triplet structure is clearly shifted with 
respect to observed frequencies for both rotational velocities. In 
the case of the 1 .70 model, the predicted m - frequency 
obtained when considering Q = 50kms ' is centred with re- 
spect to the observed ones. However the m - +\ branches are 
still far from the observations. In contrast, when considering 
the Q. - 30kms ' rotational velocity, theoretical frequencies 
tend to the observed ones for the low effective temperatures. 



Such results would thus suggest that 1 .70 pseudo-rotating 
models, with effective temperatures near the lower limit of the 
error box, and with a rotation velocity around 30 km s"' may be 
considered as representative for HD 48501. In addition, the ro- 
tational velocity confirms the spectroscopic vsin; observation. 

In order to confirm the triplet possibihty, it is worth to anal- 
yse the behaviour of the average splitting predicted by mod- 
els. For the rotational velocities considered here, the frequency 
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Fig. 10. Mean splitting (left panel) and asymmetry (right panel) of the selected {{,n) - (1,22) triplet of the l.TOM© models as 
calculated in Eqs. ^] and respectively. Solid lines represent observed values. Dotted and dashed lines represent theoretical 
values obtained from models with a rotational velocity of 30 and 50 km s ' respectively. 



structure of a rotationally split [ - I can be studied by using 
the mean splitting form 

As = ^(tr+i -fr_i) . (16) 

In Fig. [lO| (left panel) the evolution of the theoretical splitting 
{{ - \, g2i) is shown for the model giving the best frequency 
results (I.7OM0). The solid line represents the observed aver- 
age spacing, which value is A°''^ = 1.14/iHz. When consider- 
ing a rotational velocity of 30 km s ' (dotted line), the observed 
mean splitting is reproduced (A* -A°''^ ~ lO"-'- 10"^) by mod- 
els in the region of low temperatures. 

In addition, the study of the observed asymmetry may 
furnish additional information to identify the corresponding 
triplet. An equivalent form to Eq can be used to analyse 
the corresponding I - \ triplet asymmetry: 

Aa = 0-+1 -I- cr_i - 2cro , (17) 

where the subscripts -1-1,-1 and indicate the corresponding 
values of the azimuthal order m. The evolution of theoretical 
Aa values as obtained from models rotating with 50 km s"' and 
30kms"' are compared with the observed one AjJ^ = -0.16 
in Fig. [lO| (right panel). Our best predicted asymmetry is ob- 
tained for models with rotational velocities of 30 km s ' . In the 
region of low temperatures these models show A*'* approach- 
ing to the observed A*'* value in approximately 10"^ orders of 
magnitude. 

The best results for Aa and As are given for a 1.70 M© 
model with a rotational velocity of f2 = 28.3 km s"' , a radius of 
2.43 R0, and presenting an effective temperature of 6839.24 K 
(log Teff = 3.835) and log g - 3.896. Its central hydrogen frac- 
tion, Xc - 0.26, shows an evolutionary stage corresponding to 
a main sequence star with 1 .4 Gyr. 

6. Conclusions 

In the present work, the Frequency Ratio Method (FRM) pro- 
posed bv iMova et al.. (.2005i) is revisited in terms of the appli- 



cability to rotating y Doradus stars. The accuracy of the FRM 
in presence of rotation has been exhaustively examined. To do 
so, the following aspects have been considered: the validity of a 
perturbative approach to compute adiabatic oscillation frequen- 
cies of g modes in asymptotic regime; the effect of rotation on 
the observational Brunt-Vaisala frequency integral and, finally, 
the problem of multiplet-like structures. This problem concerns 
the possibility of disentangling whether the observed frequen- 
cies belong to any rotational split multiplet or to the period 
spacing expected for high-order gravity modes in asymptotic 
regime. 

It is found that reliable results may be obtained when ob- 
jects rotate with Q < 70kms ^ In such cases, the error of the 
FRM increases on order of magnitude respect to the typical er- 
rors given in paper I. 

Concerning the multiplet-like structures, all the possible 
confusing scenarios have been investigated. Provided any ad- 
ditional information on the mode degree C, simulations indi- 
cate that the FRM would be discriminating for m - modes. 
For such cases, any misinterpretation induced by the pres- 
ence of rotationally split multiplet-like structures is avoided. 
When ^ is a priori unknown, such discrimination is not en- 
sured. Nevertheless, we show here that the possibility of a con- 
fusing scenario is rather unlikely. This constitutes a very im- 
portant result, since it reinforces and extends (up to some ex- 
tent) the applicability of the FRM to slowly-moderately rotat- 
ing y Doradus stars. 

In order to test the present results with a real star the FRM is 
applied to the y Doradus star HD 48501, which can be consid- 
ered as a slow rotator (vsin/ - 29 km s ' ). The models fulfilling 
the coiTesponding («,, 3) predictions were found to be incon- 
sistent with the observational constraints for i - 1 modes, in 
particular with the observed metallicity. However, for t - 2 
modes, the FRM predicts models with masses from 1.4 to 
1 .6 Mq and radial order ranges in the range n - [38, 47] . 
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Fu rthermore, in contrast to predictions given bv lAerts et alJ 
ll2004 . when analysing the observed frequencies as belonging 
to a rotational induced ( - I triplet, solutions reproducing the 
observed average splitting and the triplet average asymmetry 
have been found to be plausible for a 1 .70 Mg model. 
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